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1. Introduction 

Building scalable multi-qubit systems is presently the main challenge towards the 
implementation of a solid-state quantum information processor [l]. High-fidelity single 
qubit gates based both on semi- [2] and super-conducting technologies are nowadays 
H ' available. In particular, for the three basic types of superconducting qubits (charge, 
flux and phase) single-qubit operations with high quality factor have been demonstrated 
in different laboratories |3l H]. Further improvement has been recently achieved with 
circuit-QED architectures [5l|6]. Multi-qubit systems instead have been proved harder 
to operate, the main limitation arising from the broadband solid state noise. The 
requirements for building an elementary quantum processor are in fact quite demanding 
on the efficiency of the protocols. This includes both severe constraints on readout and 
a sufficient isolation from fluctuations to reduce decoherence effects. 

Based on experiments on the different Josephson junction (JJ) setups, there is 
presently a general consensus on the most common noise spectral forms and on the 
main consequences on systems evolutions. Typically noise is broadband and structured. 
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i. e. the noise spectrum extends to several decades, it is non-monotonic, sometimes a 
few resonances are present. 

Noise with 1/f spectrum is common to virtually all nanodevices. Its physical 
origin varies from device to device and depends on the specific material. Different 
implementations are in fact more sensitive either to charge, flux or to critical current 
fluctuations with spectral density scaling as the inverse of the frequency. The presence 
of slow components in the environment makes the decay of the coherent signal strongly 
dependent on the experimental protocol being used IH [3 El [9]. Measurements 
protocols requiring numerous repetitions are particularly sensitive to the unstable 
device calibration due to low-frequency fluctuations. The leading effect is defocusing 
of the measured signal, analogous to inhomogeneous broadening in NMR [10]. The 
intrinsic high-frequency cut-off of 1// noise is hardly detectable, measurements typically 
extending to 100 Hz (recently charge noise up to 10 MHz has been detected in a 
SET [Hj). Incoherent energy exchanges between system and environment, leading 
to relaxation and decoherence, occur at typical operating frequencies (about 10 
GHz). Indirect measurements of noise spectrum in this frequency range often suggest 
a white or ohmic behavior [121 IE]- In addition, narrow resonances at selected 
frequencies (sometimes resonant with the nanodevice relevant energy scales) have being 
observed [TH [15] . In certain devices they originate from the circuitry [16] and may 
eventually be reduced by circuit design. More often, resonances are signatures of the 
presence of spurious fluctuators which also show up in the time resolved evolution, 
unambiguously proving the discrete nature of the noise sources [I?]. Such fluctuators 
may severely limit the reliability of nanodevices [9l [18] . 

Explanation of this rich physics is beyond phenomenological theories describing 
the environment as a set of harmonic oscillators. On the other side, an accurate 
characterization of the noise sources might be a priory inefficient, since a huge number of 
parameters would be required for a microscopic description. Therefore one may follow a 
different route, consisting in classifying noise sources on the basis of their effects instead 
of on their nature and to understand, case by case, which is the efficient description 
of the environment. The required information may in fact depend on the specific 
protocol, especially when the environment is long-time correlated. This program is 
meaningful for quantum information where relevant time scales are much smaller than 
the decoherence time. This means that in the generic favorable situation coupling 
with the environment has simple effects on the system dynamics. As a difference 
with condensed matter physics where long time behavior is emphasized and interesting 
problems involve entanglement of a many-body system. 

Here we illustrate a road-map to treat broadband noise which allows to obtain 
reasonable approximations by systematically including only the relevant information on 
the environment, out of the huge parametrization needed to specify it. The multi-stage 
approach for the different classes of broadband noise has been originally introduced for 
single qubit gates [9]. The obtained predictions for the decay of the coherent signal 
are in agreement with observations in various JJ implementations and in different 
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protocols [SI in]. We mention the observed decay of Ramsey fringes in charge-phase 
qubits jl] and recent results on flux/phase qubits |T9]. Here we extend the procedure 
to complex solid-state architectures. As an illustrative case, we perform a systematic 
analysis of the effects and interplay of low and high frequency noise components in a 
two-qubit gate in a fixed coupling scheme. Such a systematic analysis points out that 
efficient operations in the solid state require an accurate preliminary characterization of 
the noise spectral characteristics and tuning appropriately the device working point. 

The paper is organized as follows: in Section 2 we introduce the general framework 
of solid-state multi qubit gates and illustrate the characteristics of the most common 
spectral forms. In Section 3 we propose a classification of the noise sources and in 
Section 4 we present a multi-stage approach to deal with broadband noise. Section 5 
is dedicated to adiabatic noise. In Section 6 the multistage approach is applied to a 
universal two-qubit gate. Section 7 summarizes our main findings. 

2. Multi-qubit systems and noise 

A multi-qubit gate can be modeled by the following Hamiltonian [h = 1) 

nG{t) = J2n^^{t) + Y.'^.At) (1) 

where for each qubit, labeled by z = l,...,n, T-L^Qit) = —\ ^i{t) ■ (t^*^ includes 
set of parameters intrinsic to the device and time dependent classical control fields. 
Manipulation of tunable fields allows control of the dynamics and design of arbitrary 
single qubit operations via unitary transformations. Interactions among qubits and 
the needed additional parametrization to describe the multi-qubit gate are included in 
'^ij'Hijit). Depending on the design, the control Hamiltonian may span a reduced 
subspace of the qubit Liouville space. In other words, there is a limited number of ports 
available for control, for tuning, for state preparation and measurement. Noise is also 
coupled via these ports. It is rather usual that only one of the control fields allows the 
required fast addressing of each qubit. Therefore we focus on a model where both the 
control fields, A{t), and the environment are coupled to a single qubit operator, say az\ 

n§ {t)= - ^ cos a« - ^ sin 6, a« - i a« , (2) 
nit) = HcW + i + (3) 

i 

where the polar angles 6i define qubit-i working point. The environment Hamiltonian 
is Hji and Xi is a collective environment variable acting on qubit i. Model implies 
a projection of the device Hamiltonian onto the subspace spanned by the two lowest 
energy eigenstates for each qubit. This description is valid provided that manipulations 
with the control fields do not induce leakage to higher energy states of the device |20j . 

We leave unspecified the nature of the noise sources described by H/j, relevant cases 
being either discrete or Gaussian fluctuations. In addition, the environment Hamiltonian 
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may include correlations among noise sources affecting each qubit. In the spirit of the 
present analysis, we assume that the only information at our disposal is the power 
spectrum of Xi fluctuations 

SxA^) = dt e^' 1^ (x,(t)x,(o) + x,(o)x,(t)) - (x,)2| , (4) 

where (...) denotes the equilibrium average with respect to 'Kr. We remark that the 
effect of the environment on the system dynamics is completely characterized by Sx^ (i^) 
only if the environment is composed of harmonic oscillators or if it is weakly coupled 
to the system and short time correlated. Of course this is not the case if low-energy 
excitations determine memory effects. This is the typical situation in the solid state 
where in general additional statistical information on the environment is required. 

We assume that each noise component Xj has broadband spectrum SxX^) — ^5 
uj G {^m-,lM] followed by a white or ohmic flank at frequencies uj ^ 7^. Low- 
and high-frequency cutoffs depend on the specific setup. Impurities of various origin 
responsible for random telegraph fluctuations contribute Lorentzian peaks to the power 

2 

spectrum 5'rtn(i^) = ■^'^Xufi ■ Such a spectrum originates from classical fluctuations 
X — )■ X{t\ where X{f) randomly switches between two values {0,fo} with rate 
7. An ensemble of bistable impurities with a distribution of switching rates 
oc 1/7, 7 G [7m)7Af] and average coupling strength gives rise to l//-spectrum, 
S^l^{uj) ~ [(7r/4) Xfcjt>2/ ln(7jv//7m)] f^""*^ [21]. Selected Lorentzian peaks may be 
visible in the spectrum if individual impurities are strongly coupled, i. e. when 
fo/7 ^ 1- Instead, damped coherent fluctuators at frequencies Wq; in the simplest 
cases contribute to the total spectrum with additional peaks, Scf{oS) = ^a^(2-a;o)2 • 
Spectroscopic evidence of coherent impurities of frequency ujq close to the qubit Larmor 



frequency has been reported in [Ml [22]. Remarkably, the possibility to exploit spurious 
quantum two-level systems as qubits [23] or for quantum memory operations has been 
recently demonstrated [21]. In these cases a quantum description of the impurity is 
required 



3. Three classes of noise 



The above description illuminates that in the solid state we have to deal with broadband 
and structured noise. In other words, the noise spectrum extends to several decades, it 
is non-monotonic, sometimes a few resonances are present. The various noise sources 
responsible for the above phenomenology have a qualitative different influence on the 
system evolution. This naturally leads to a classification of the noise sources according 
to the effects produced rather than to their specific nature. 

The effects of high-amplitude noise at low frequencies, like 1// noise, vary from 
protocol to protocol. This feature is typical of non-Markovian baths. Quantum 
operations necessarily require repetitions of single detections, each leading to Boolean 
answer. Therefore, even "single shot" measurements result from numerous repetitions of 
single runs in an overall process which may last minutes. In the presence of low frequency 
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fluctuations this leads to unstable device calibration and random clock frequencies in 
the various repetitions. As a result, a de-focused signal is observed, a phenomenon 
analogous to inhomogeneous broadening in NMR. Re-focusing protocols, like echo 
or some dynamical decoupling schemes, allow partial recovery of the signal [TJ [13]. 
Since the environment is long-time correlated, statistical information beyond the power 
spectrum may be required to describe its effects. This is the case for instance in echo 
protocols [8l|18]. Environments with long-time memory belong to the class of adiabatic 
noise, for which the Born-Oppenheimer approximation is applicable. We classify this 
part of the noise spectrum as " adiabatic noise" . 

Noise at higher frequency, around the system typical scales (e.g. single qubit Larmor 
frequencies), results in incoherent energy exchanges between the quantum device and 
the environment. In particular, high-frequency noise is responsible for spontaneous 
decay. For systems relevant for quantum information, the system-bath coupling is rather 
small, in addition high-frequency noise is short-time correlated. Therefore, in simplest 
cases, effects can be described by a Born-Markov master equation [27]. For single qubit 
gates it leads to the relaxation and secular dephasing times, Ti = sin^ ^S'jjf (n)/2 and 
T2 = 2Ti [To] . We classify this part of the noise spectrum as " quantum noise" . 

Finally, resonances in the spectrum unveil the presence of discrete noise sources 
which severely effect the system performances, in particular reliability of devices. This is 
the case when classical impurities are slow enough to induce a visible bistable instability 
in the system intrinsic frequency. For instance, single qubit gates in the presence of 
random telegraph noise (switching rate 7) may display two effective frequencies, Q 
and Q', depending on the impurity state. Their visibility is measured by the ratio 
g = [Q' — Q)/'y. Beatings can be observed in the "strong coupling" regime g > 1 [28]. 
Quantum impurities may also entangle with the device. This additionally leads to a 
variety of features, like peculiar temperature dependencies of decay rates [25]. Under 
these conditions, knowledge of the power spectrum is absolutely insufficient and in order 
to describe these effects the relevant system Hilbert space has to be enlarged to include 
the responsible environmental degrees of freedom. Effects in general depend on the 
specific protocol and require a microscopic model of the fluctuators. We classify this 
part of the noise spectrum as "strongly-coupled noise". 

Each noise class requires a specific approximation scheme, which is not appropriate 
for the other classes. The overall effect results from the interplay of the three classes of 
noise. In the following Section we will illustrate a multi-scale theory to deal with solid 
state broadband noise. 

4. Multi-scale theory for broadband noise 

We are interested to a reduced description of the n-qubit system, expressed by the 
reduced density matrix p"'{t). It is formally obtained by tracing out environmental 
degrees of freedom from the total density matrix W^'^'^'~^(t), which depends on quantum 
(Q), adiabatic (A) and strongly coupled (SC) variables. The elimination procedure can 
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be conveniently performed by separating in the interaction Hamiltonian, ^jCri*'' ® Xi, 
various noise classes, e.g. by formally writing 

ai'^ ® X, = ® X^ + af> ® Xf + ® Xf ^ . (5) 

Adiabatic noise is typically correlated on a time scale much longer than the inverse of 
the natural frequencies then application of the Born-Oppenheimer approximation is 
equivalent to replace Xf with a classical stochastic field Xf{t). This approach is valid 
when the contribution of adiabatic noise to spontaneous decay is negligible, a necessary 
condition being t <^ oc S'^(fij)^^. This condition is usually satisfied at short enough 
times, since Sx{(-o) is substantially different from zero only at frequencies u <t^Qi. 

This fact already suggests a route to trace-out different noise classes in the 
appropriate order. The total density matrix parametrically depends on the specific 
realization of the slow random drives X{t) = {Xf{t)} and may be written as 
W'^''^'^'-^(t) = W'^'^'-^(t\X{t)). The first step is to trace out quantum noise. In the 
simplest cases this requires solving a master equation. In a second stage, the average 
over all the realizations of the stochastic processes, X{t), is performed. This leads 
to a reduced density matrix for the n-qubit system plus the strongly coupled degrees 
of freedom. These have to be traced out in a final stage by solving the Heisenberg 
equations of motion, or by approaches suitable to the specific microscopic Hamiltonian 
or interaction. For instance, the dynamics may be solved exactly for some special 
quantum impurity models at pure dephasing, 9i = 0, when impurities are longitudinally 
coupled to each qubit [HI EH] . The ordered multi-stage elimination procedure can be 
formally written as 

p-{t)=Trsc^j P[X(t)]P[X(t)]TrQ[iyQ'^^(t|X(t))]| . (6) 

In the following Section we concentrate on the elimination of adiabatic noise. The 
ordered procedure will be illustrated in Section 6 for a two-qubit gate. 



5. Adiabatic noise 



In general, the adiabatic approximation holds true for times short enough to fulfill the 
necessary condition, t <^ T^. In the peculiar, pure dephasing regime, 6i = 0, relaxation 
processes are forbidden and the adiabatic approximation is exact for any SxX^)- 
addition, the adiabatic scheme can be applied also in the presence of correlations between 
processes Xi(t) and Xj(t) [29j. For the sake of clarity here we consider adiabatic noise 
affecting independently each qubit. Moreover, we exclude time-dependent drives in A{t) 
and in 'Hij{t). The procedure can be straightforwardly extended for instance to Rabi 
oscillations and other ac-drives [30ij . 

Suppose we are able to diagonalize, for each realization X{t) of the stochastic 
processes, the system Hamiltonian, 

^(^) = -Y.\^^- ^^'^ + E + ^ E ® ^^(t) , (7) 

i i,j i 
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and denote \m{Xt)) an instantaneous eigenstate of ([7]) with eigenvalue Emit). If 
the system is prepared in a pure state, | ^Z' ) , in the adiabatic approximation, each 
component | m(Xo) ) ( m(Xo) \ ipQ) evolves in time according to 

|m(X<))e'*™(*)(m(Xo)|V^O) , 

where ^^{t) = —i j^dsEjn{s) is the dynamic phase. The stray fields X(t) affect the 
time evolution of the gate fidelity via modifications of the phases of the 

eigenstates | m(Xf) ) . In addition, | m(Xo) ) reflects imperfect preparation due to the 
random field. This occurs, for instance, when the system is prepared by a reset to a 
pure state vector | ^/^ ) followed by an initialization pulse whose effect depends also on 
the stray field at t = 0. The evolution of the system conditional density matrix can be 
presented in a compact form 

p"(t|X(s)) = 5^i?^,[Xo,Xi] e-^/o^^f^™.(^) , (8) 

mp 

where the instantaneous splittings appearing in the phase factor are VLmp{t) = 
— Qpm(t) = Emit) — Ep(t) and we have introduced the operator 

Rmp[Xo,Xt]= \m{Xt)) (m(Xo)|p(0)b(Xo)) {p{Xt)\ (9) 

which contains information about the preparation and the eigenstates errors. Finally, 
the average of p^{t\X{s)) over the stochastic process yields the system density matrix, 
which can be presented path integral 

p-{t) = I V[X{s)] P[X{s)] p-{t\X{s)) . (10) 

Here P[X(s)] contains information both on the stochastic processes and on details of 
the specific protocol. It is convenient to split it as follows 

P[X{s)] = F[X{s)] p[X{s)] , 

where p[X(s)] is the probability of the realization X{s) . The filter function F[X(s)] 
describes the specific operation. For most of present day experiments on solid- 
state qubits F[X{s)] = 1. For an open-loop feedback protocol, which allows initial 
control of some collective variable of the environment, say Xq = 0, we should put 
F[X{s)] oc 6{Xo) = n,5(X,o). 

A critical issue is the identification of p[X(s)] for the specific noise sources, as those 
displaying 1// power spectrum. If we sample the stochastic process at times t^ = kAt, 
with At = t/m and = 0, . . . , m, we can identify 

p[X{s)]= lim p„+i(X„,t;...;Xi,ti;Xo,0), (11) 

m— >oo 

where Pm+i{-) is a m + 1 joint probability and we have used the shorthand notation 
Xk = Xt^.. In the following we will propose a systematic method to select only the 
relevant statistical information on the stochastic process out of the full characterization 
included in p[X(s)]. 
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We would like to remark, at this point, that in the adiabatic treatment the word 
decoherence is perhaps abused. Decoherence is ultimately due to entanglement of the 
qubit to a quantum environment [31l |32] , whereas classical adiabatic noise produces only 
de-focusing. However, unless the signal can be totally re-focused, and this is impossible 
in practice, the behavior of p{t) is the same as for true decoherence [33]. In other words, 
although proper and improper mixed states at the fundamental level are well distinct 
concepts [31], they are not distinct in the density matrix description. We also observe 
that present day "applied" research on solid-state coherent nanodevices for quantum 
information focuses on the " short-time" dynamics, since a signal which is almost decayed 
is useless. In this context, methods as the adiabatic approximation are valuable even if 
they are not accurate at long time scales and are not valid down to zero temperature. 

5.1. Longitudinal approximation 

The longitudinal approximation consists in neglecting modifications of the eigenstates 
I m{Xt) ) and preparation effects. Without loss of generality, we may assume vanishing 
average of the stochastic processes after the preparation pulse, Xq. The longitudinal 
approximation amounts to put \m{Xt)) = \m{Xo)) = |m) where \m) is an 
eigenstate of the system Hamiltonian "Hg- Then Rmp[Xo, Xt] ~ l'^)Pmp(O) {p\ is a 
projected element of the initial density matrix and Eq.l fTOj) simplifies to 

Pipit) = p-^^{0) I V[X{s)] P[X{s)] e-/o^^^-W . (12) 

The significance of the longitudinal approximation is easily illustrated for a single qubit 
gate, when the system Hamiltonian reduces to 

nit) = COS 0a, - I smea,-^a,®X{t). (13) 

In this case, retaining variations of the splitting Q amounts to consider only fluctuations 
of the length of the vector Q. This would be the only effect if the noise would act 
longitudinally, 6 = 0. Variations of the eigenstates | m{Xt) ) originate instead from 
"transverse" variations of Q . 

The longitudinal assumption has been performed in Ref. ^35j to discuss the effect 
of Gaussian adiabatic environments. The present approach automatically provides 
constraints on its validity and shows that whereas errors due to transverse fluctuations 
are weakly dependent on time, phase errors accumulate. Therefore, transverse noise in 
the adiabatic approximation has possibly some effect only at very short times, but the 
phase damping channel eventually prevails. These considerations strongly depend on 
the amplitude of the noise. We checked analytically and with simulations that they hold 
true for realistic figures of noise as those measured in experiments by Zorin et al. [36] . 

In the longitudinal approximation, diagonal elements of the reduced density matrix 
in the eigenstate basis, | m ) = | ± ) , do not decay. Instead, the decay of the off- 
diagonal elements results from p it) = p (0)exp[— ^^(t)], where the complex average 
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phase is given by 



^{t) = -nt + i\n V[X{s)]P[X{s)]e'-^od^^[x^'^K 



(14) 



We notice that the longitudinal approximation may be exact for certain protocols, which 
therefore are not affected by adiabatic transverse noise. For instance, this is the case 
if the system is described by ( IT3l) and it is prepared in a state p(0) = ^(1 ± ay). If we 
measure ay, since {m{Xt) \ ay \ p{Xt) ) do not depend on Xt, the result is not affected 
by transverse fluctuations and {ay{t)) = {ay{0)) cos[7?.e$(t)] exp{Xm$(t)}. This the 
case of the decaying oscillations pattern measured with Ramsey interference [1], if the 
effect of imperfect 7r/2 pulses is negligible. 

5.1.1. Static path approximation A standard approximation of the path integrals 
fll0|12l) consists in neglecting the time dependence in the path, X{s) = Xq and taking 
F[X] = 1. In this Static Path Approximation (SPA) the problem reduces to ordinary 
integrations with Pi{Xq, 0) = p{Xq). In the single qubit case, for instance, Eq. flT^ gives 
average phase shift 



where fi(Xo) = ^y[^l sin 6 + Xq]^ + (fi cos 6')2. Eq.(^ describes the effect of a 
distribution of stray energy shifts flrnpi^) — ^mp(O) and corresponds to the rigid lattice 
breadth contribution to inhomogeneous broadening [10]. In experiments with solid state 
devices this approximation describes the measurement procedure consisting in signal 
acquisition and averaging over a large number N of repetitions of the protocol, for 
an overall time tm (which may also be minutes in actual experiments). Due to slow 
fluctuations of the solid state environment calibration, the initial value, f2sin^ + Xq, 
fluctuates during the repetitions blurring the average signal, independently on the 
measurement being single-shot or not. 

The probability p{Xq) describes the distribution of the random variable obtained 
by sampling the stochastic process X{t) at the initial time of each repetition, i. e. at 
times tk = ktm/N, k = 0, N — 1. If Xq results from many independent random variables 
of a multimode environment, the central limit theorem applies and p{Xq) is a Gaussian 
distribution with standard deviation a 



with integration limits 1/tm, and the high-frequency cut-off of the 1// spectrum, •jm- In 
the SPA the distribution standard deviation, a, is the only adiabatic noise characteristic 
parameter. If the equilibrium average of the stochastic process vanishes, Eq. (1T5|) reduces 
to 




(15) 





(16) 
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A convenient approximation is obtained by expanding Q{Xq) to second order in Xq, 
which leads to |9] 

The short-times decay of coherent oscillations qualitatively depends on the working 
point. In fact, the suppression of the signal, exp[Xm$(t)], turns from a 
exp(— i(cos6'(Tt)^) behavior at 6* to a power law, [l + (sin6'^(j^t/fi)^]~^/^, at 6' pa 7r/2. 
In these limits Eq. (fT7|l reproduces known results for Gaussian 1// environments. In 
particular, at ^ = we obtain the short-times limit, t <^ 1/7a/, of the exact result of 
Ref . [32] . At = 7r/2 the short and intermediate times result of Ref. [35] is reproduced. 
The fact that results of a diagrammatic approach with a quantum environment, as those 
of Ref. [35] , can be reproduced and generalized already at the simple SPA level makes 
the semi-classical approach quite promising. It shows that, at least for not too long times 
(but surely longer than times of interest for quantum state processing), the quantum 
nature of the environment may not be relevant for the class of problems which can be 
treated in the Born-Oppenheimer approximation. Notice also that the SPA itself has 
surely a wide validity since it does not require information about the dynamics of the 
noise sources, provided they are slow [37] . 



5.1.2. Beyond SPA: first correction Going beyond the SPA amounts to sample more 
accurately the stochastic process X{s) in ( fTTl) . The first correction to the SPA 



t 



s. 



is obtained by parametrizing the random process as follows X{s) = Xq + 
Inserting this expression in (fTOj) and approximating p[X(s)] in (fTTj) with the conditional 
probability p[Xt, t; Xq, 0], we obtain 



Pit) 



dXt dXoP2{Xt, t; Xo, 0) p(t|Xo + ^^% ) 



The joint probability depends on the statistics of the noise sources. Therefore the first 
correction to the SPA distinguishes discrete and Gaussian processes. Considering again 
a single qubit, the average phase ( IT^ at the working point 6* = 7r/2 for generic Gaussian 
noise becomes 



2$(t) 



iQt + - In 
2 



a^[l-n{t)]t 



I, 
+ - In 
2 



a'^Tx{t)t 
3fi . 



where Ti{t) = ^ {du/ir) S{u){l — e *'^*) is a transition probability, depending on the 



stochastic process. For Ornstein-Uhlenbeck processes it reduces to the result of Ref. 
The first correction suggests that the SPA, in principle valid for t < I/'Jm, niay have 
a broader validity. This is illustrated in Fig.(IT]) where the adiabatic approximation 
(numerical evaluation of Eq. ffMj) ) is compared with the exact (numerical) dynamics 
in the presence of I// noise and with the analytic forms resulting from the SPA and 
its first correction. For I// noise due to a set of bistable impurities the SPA valid 
also for t 3> 1/7m, if 1m ^ ^- Of course the adiabatic approximation is tenable if 

t<Tf = 2/s^{n). 




Figure 1. Imaginary part of tlie colierence (t) at 8 — 7r/2: Numerical simulations 

of an adiabatic 1// environment, S'^/^{uj) ~ [{tt / 4) Ni,fv'^ / \\i{'^M /im)] , 
correspond to the thin black line. Noise is produced by = 250 bistable fluctuators 
per decade, with = 10^ rad/s < 7i < 7m = 10^ rad/s < Q. ^ 10^° 

rad/s. The coupling v = 0.02 is appropriate to charge devices, and corresponds to 
S{ijj) — 16ttAE'^/uj with A ^ 10~^ [35]. The adiabatic approximation fully accounts 
for dephasing (numerical evaluation of Eq. p^ . red dot-dashed line). The Static Path 
Approximation (SPA) Eg. dTTl) (blue solid line) and the first correction (blue dashed 
line) account for the initial suppression, and it is valid also for times i S> 1 /7m • 



By sampling more accurately the adiabatic process X{t) it is possible to selectively 
include the statistical information needed for the specific measurement process. For 
instance, echo protocols are able to partly re-focus the signal, in other words de-focusing 
described by the SPA is almost canceled by the echo pulse sequence. The decay of the 
echo signal is due to the un-canceled dynamics of the low-frequency fluctuations and 
to quantum noise. The leading effect of adiabatic noise can be estimated by a proper 
parametrization of X{t), similar to the one considered in the present paragraph [30]. 
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6. Two-qubit universal gate: multi-stage approach 

In the present Section we apply the multi-stage approach to a universal two-qubit gate 
based on a fixed coupling scheme. Capacitive and inductive fixed couplings |39] have 
been used to demonstrate two-qubit logic gates in different J J implementations [40] . 
Entanglement is generated by tuning single-qubit energy spacing to achieve mutual 
resonance. To this end, during the gate operation at least one qubit has to be 
moved away from the working point of minimal sensitivity to parameters variations, 
the "optimal point". This has so far represented the main drawback of the fixed- 
coupling scheme for Josephson implementations, with the exception of phase qubits |41j . 
More recent proposals have attempted to solve this problem by introducing tunable 
coupling schemes |12]. Most of them rely on additional circuit elements and gain their 
tunability from ac-driving [13] or from "adiabatic" couplers Some of these coupling 
schemes have been tested in experiments and are potentially scalable [15]. None of these 
implementations is however totally immune from imperfections. In general, introducing 
additional on-chip circuits elements opens new ports to noise. The possibility to employ 
"minimal" fixed coupling schemes has been recently reconsidered in Ref. [46], pointing 
out the possibility to single out "optimal coupling" conditions which ensure reasonable 
protection from 1// fiuctuations. 

In order to implement a y/i — SWAP gate in a fixed coupling scheme we need two 
resonant qubits with a transverse coupling, as described by the Hamiltonian 

no = ^i^) I(^) - § I« -i^) + I ai^^ ® ai^) (19) 

where Uc is the coupling strength, and ai* ■* the pseudo-spin operators whose eigenstates 
|±) (eigenvalues ±1) are the computational states of qubit i. Eigenvalues and 
eigenvectors of T-Lq are given by 

uo= -Q^l + gy4 , |0) 
cui = - uj2 , |1) 

a;2=Wc/2 , |2) 

uj^ = ny/l + gy4 , |3) 

where simp = (7/(2^/1 + g'^/4:), cos(p = — 1/a/1 + with g = ooc/^ and we have 
used the shorthand notation = ® |z/)2 fJ'ji^ & {+)"}■ The level structure of 
the coupled system is schematically illustrated in Fig. ([2]). As a result of the diagonal 
block structure of the Hamiltonian f|T9|) in the computational space, the two-qubit 
Hilbert space is factorized in two subspaces spanned by pairs of eigenvectors. A system 
prepared in | H — ), freely evolving for a time tE = vr/2(x'c, yields the entangled state 
(I H — ) — "^1 — corresponding to a a/I — SWAP operation. The pair of states 
|1) and |2) span the subspace where the \/i — SWAP gate is realized, which we name 
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Figure 2. (a) Level structure of the uncoupled resonant qubits; (b) levels in the two- 
qubit Hilbert space and logic basis of product states, (c) When the coupling is turned 
on |H — ) and | — h) mix and an energy splitting Wc ^ develops between the eigenstates 

{|2),|1)}, spanning the SWAP subspace. Product states | ) and | + +) weakly 

mix and split, with 2^l^/l + g^/4:. The eigenstates {|0),|3)} span the Z subspace. 
Transverse noise ()24[) only mix eigenstates belonging to the two subspaces. Dissipative 
transitions between subspaces are indicated by wavy lines. At typical temperatures, 
ksT <C il, thermal excitation processes between subspaces can be neglected. 



SWAP-subspace. The subspace spanned by the pair of states |0) and |3) is termed Z- 
subspace. Eq. (fT9|) models, for instance, capacitive coupled charge qubits or inductively 
coupled flux qubits [40j. Typically, these system can be driven by pulses actiiig along 
the transverse direction, ai^\ The same port introduces noise into the system Here 
we consider noise sources acting transversely with respect to each qubit, i. e. the 
interaction Hamiltonian takes the form 

Ui = a(^) ® I2 Xi - ^ Ii ® ) X2. (24) 

Each noise component Xi has broadband spectrum SxX^) = ^, u E {7m,7Af} followed 
by a white flank at frequencies u > 'Jm- Correlated noise sources acting on both qubits 
have been addressed in [29] . 

If the system is initialized in the SWAP subspace, for instance in the state | H ), 

bitwise readout gives the qubit 1 switching probability P^^\t), i.e. the probability that 
it will pass to the state |— ), and the probability P^^)(t) of finding the qubit 2 in the 
initial state |— ). Cyclic anti-correlation of the probabilities signals the formation of the 
entangled state during the a/I — SWAP operation. In terms of the two qubit reduced 
density matrix in the eigenstate basis the switching probabilities read 

P^'\t) = (-|Tr2p(t)|-) = ^ [pn{t)+p22{t)]+Poo{t) + [psaW - PooW] sin^ | 

+ Re[pi2(t)] + Re[po3(t)] sinyp (25) 

I In principle noise can also couple longitudinally, i. e. via cri*'' IMj- This is the case of the charge-phase 
two-port design of Ref. [1] . 
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20 40 60 80 100 

t (ns) 

Figure 3. Interplay of adiabatic and quantum noise in the switching probability of 
qubit 1. Adiabatic noise in the SPA is parametrized by cr/fJ = 0.08. The effect of 
low-frequency noise is shown in red, the effect of white noise, Sxii^^) ~ 8 x lO^s""'^, 
in gray. Interplay of adiabatic and quantum noise is shown by the black line. The 
interaction strength is uJc/^ = 0.01. 

P^^\t) = (-lTrip(t)l-) = ^ [pii(t) +P22(t)] +Poo(t) + [P33(t) -poo(t)]sin2| 

- Re[pi2(t)] + Re[po3(t)] simp. (26) 

The matrix elements entering the above probabihties can be evaluated in the multi-stage 
approach. 

6.1. Multi-stage approach 

We split the interaction as in Eq.(|5]), Xi = Xf^ + Xi{t). Low- frequency noise is treated 
in the adiabatic and longitudinal approximation. In addition we limit the analysis to 
the SPA. We denote with uJi{Xi,X2) the eigenvalues of "Hq + "^i with Xi — )• {Xi}. 

First stage: elimination of quantum noise Quantum noise is traced out by solving 
the Born-Markov master equation for the reduced density matrix [49] with eigenvalues 
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parametrically dependent on the random fields {Xj}. In the system eigenstate basis 
and performing the secular approximation (to be self-consistently checked) it takes the 
standard form |27j : 

Pijit) = - i^^ij + Tij) Pij{t) , (28) 

where Uij = cOi{Xi, X2) — cOj{Xi, X2). The rates Tim, Tjj depend on the real parts of the 
lesser and greater Green's functions, describing emission (absorption) rates to (from) 
the quantum reservoirs |17]. In addition, for white quantum noise energy shifts are 
vanishing and do not appear in Eq. f l28|) . Because of the symmetry of (|T9|) . dissipative 
transitions inside each - SWAP or Z - subspace are forbidden. The allowed inelastic 
energy exchange processes are evidenced in Fig. ([2]). The only independent emission 
rates are Fio = r32, = T^i. They read 

Tio = 1(1 + sin(^) [Cxi(wio) + Cx^A^io)] ^^29) 
= I (1 - sinv9) [Cxi(w2o) + C'x2(w2o)] , 

2 

where the absorption rates, CxX^) = i^e^p (-uj/kBT) ' related to the spectrum of 
quantum noise, S^.{uj). Emission rates have the same form with CxX^^im) replaced by 
C Xi{—^im) ■ For the considered initial condition, | H — ), the only non vanishing elements 
of the reduced density matrix in the eigenbasis of 1-Lq are the populations and the SWAP 
coherence, P2i(t). For independent quantum noise sources acting on the two qubits, the 
SWAP decay rate reads = ^ [rio + roi + r2o + ro2] = | [11+12], where = Fio + Foi, 
i = 1, 2, are relaxation rates of the SWAP levels. From Eq. (j28|) we get 



P2i(t) = -i e-2i(Xi,X2)t g-ri2* _ (30) 

The secular approximation is valid provided that U21 ~ ^ F12. This condition is 
fulfilled, for instance, for white noise levels extrapolated from single charge-phase qubit 
experiments, as it can be evinced from the slow decay due to quantum noise reported 
in Fig. ([3]), gray line. 

Equations (1271) for the populations do not decouple even in the secular limit. 
General solutions are quite lengthy. Here we report the approximate forms in the small 
temperature hmit with respect to the uncoupled qubits splittings, ksT ^ Q. In this 
case level 3 remains unpopulated and 

Poo(t)~l-^[e-^^''* + e-^^°*] , Pii(t)~^e-ri«* , p22(t) ~ ^ e-^^"* . (31) 

Second stage: elimination of adiabatic noise We now consider the effect of low- 
frequencies in the adiabatic, longitudinal and Static Path approximations. Populations 
are unaffected by adiabatic noise, whereas the SWAP coherence, puit), has to be 



§ We are here disregarding effects due to the instantaneous eigenstates of Hq + Hi with Xi — > {Xi} 
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Figure 4. Effect of adiabatic noise under optimal tuning: switching probability of 
qubit 1 in the presence of low-frequency noise parametrized hy a/fl — 0.08 for "optimal 
coupling" LUc/^ = 0.08 (orange) and non-optimal tuning ujc/^ — 0.01 (red). 



averaged over Xj. Here we disregard the negligible Xi dependence of the rates Tij 
via Uij. Under these simplifying assumption, we are left with the following average 

1 



jW2l(Xi,X2)t 



(32) 



exp[—j-^]. The SWAP splitting can be estimated by treating in 



perturbation theory "Hi, with Xi — )■ {Xi}, with respect to "Hq- This leads to 

UJ2l{Xi,X2) ^ Uc 



+ 



ixf-xi 



2\2 



The average in Eq. fl32l) can be evaluated in analytic form and gives [46.J 

1 



-ri2t 



f2 2iuJc 



ei-'^t+Ht) KQ[h{t)] 



(33) 



(34) 



2^2 V vrt 

where h{t) = iUc/tiVt^ /a"^ + iujctY / {AVL^), and KQ[h] is the K-Bessel function of 
order zero We considered equal standard deviations for noise sources acting 



Broadband noise decoherence in solid-state complex architectures 



17 



on both qubits, CTj = a. Inserting (IHT]) and in (^5^ and (^E^ we obtain the 

switching probabihties in the multi-stage approach. Out of phase oscillations signals 
two-qubit states anti-correlations and follows from P^^\t) = P(t) ± Re[pi2(t)], with 
P{t) = -i cosy? [pii(t) + p22{t)] + cos2 (f ) ^ -i cosy? [e'^io* + e-^2o<] + ^2/4^ 

The efficiency of the gate results from the interplay of quantum and adiabatic noise. 
High-frequency noise levels expected from single qubit experiments [13] weakly affect 
the switching probability, whose decay is mainly due to low- frequency noise, Fig. (|3]). 
Remarkably, considerable recovery of short-times oscillation amplitude may be achieved 
by an optimal choice of the coupling strength, Uc ^ cr [3Bj . This regime is illustrated in 
Fig.©. 

7. Conclusions 

In this article we presented a road-map to treat broadband noise typical of solid 
state nanodevices. The introduced multi-stage approach allows to obtain reasonable 
approximations by systematically including only the relevant information on the complex 
environment, out of the huge parametrization which would be required for a microscopic 
description. Since the environment is in general long-time correlated, the required 
information depends on the specific protocol. 

The predictions obtained with the present approach are in agreement with 
observations in various single qubit J J implementations and in different protocols [HI [9] . 
We extended the procedure to deal with complex solid-state architectures. This is 
a required step in order to predict efficiency and possibly appropriately design of 
nanodevices for quantum information processing. Both because of the complexity of 
architectures and of the unavoidable broadband nature of solid state noise, theoretical 
tools allowing systematic and controlled approximations are particularly valuable. 

As an illustrative case, we performed a simplified analysis of the effects and interplay 
of low and high frequency noise components in a two-qubit gate in a fixed coupling 
scheme. Our results points out that efficient operations in the solid state require an 
accurate preliminary characterization of the noise spectral characteristics and tuning 
appropriately the device working point. 
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